Non-linear thermoelectrics of molecular junctions with vibrational coupling 
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We present a detailed study of the non-linear thermoelectric properties of a molecular junction, 
represented by a dissipative Anderson-Holstein model. A single orbital level with strong Coulomb 
interaction is coupled to a localized vibrational mode and we account for both electron and phonon 
exchange with both electrodes, investigating how these contribute to the heat and charge transport. 
We calculate the efficiency and power output of the device operated as a heat to electric power 
converter and identify the optimal operating conditions, which are found to be qualitatively changed 
by the presence of the vibrational mode. Based on this study of a generic model system, we discuss 
the desirable properties of molecular junctions for thermoelectric applications. 

PACS numbers: 84.60.Rb 85.65.+h, 85.35.Gv, 85.35.-p, 



I. INTRODUCTION 

The field of single-molecule electronics has been ex- 
panding rapidly during recent years, as techniques to 
electrically contact and control single molecules in a 
transport junction have improved [l|-[5|. By studying 
the electric current through the molecule as function of 
the applied voltage-bias, spectroscopic information can 
be extracted [2]. In setups with a gate-electrode, which 
can be used to control the electrostatic potential on the 
molecule, a detailed spectroscopy can be performed 
[||. By applying a temperature-bias and measuring the 
induced electric current or voltage, additional informa- 
tion can be extracted, such as the type of carriers (holes 
/ electrons) dominating transport pL This emerging 
field of molecular thermoelectrics j7Hllj is also interest- 
ing for applications. Molecules have been predicted to 
be particularly efficient for conversion of heat into elec- 
tric energy [9|, [l0| (or analogously for cooling, using elec- 
tric energy to pump heat), the reason being their very 
sharp electronic resonances when weakly coupled to elec- 
trodes [Tcj. This is similar to the large thermoelectric 
efficiency of e.g., semi-conducting nanowires with highly 
peaked densities of states [12| . 

Most theoretical works on meso- and nano-scale ther- 
moelectrics have focused on the linear, equilibrium 
regime, where one operates close to the small voltage 
V = —SAT which exactly cancels the current induced 
by the small thermal bias AT. Here the thermopower (or 
Seebeck coefficient) S = Gt/G is the decisive quantity, 
where (Gt) G is the (thermal) conductance. A large effi- 
ciency r] of the device operated as a heat to electric energy 
converter is then related to a large dimensionless thermo- 
electric figure of merit ZT = GS 2 T / where T is the op- 
erating temperature and n the thermal conductance. In 



bulk systems, ZT is normally limited by the Wiedemann- 
Franz law, stating that k/GT is a system independent 
constant. However, the Wiedemann- Franz law is a result 
of Fermi-liquid theory and breaks down in mesoscopic 
and nanoscopic systems, e.g., due to large Coulomb inter- 
action, as has been demonstrated for quantum dots [l3| 
and metallic islands [3], allowing much larger values of 
ZT to be reached. As ZT — )• oo, the efficiency approaches 
the ideal Carnot value rj rjc = 1 - T/(T + AT) pj|. 

However, in the linear regime, AT <C T, the effi- 
ciency stays low even if ZT can be made very large: 
r] ~ AT/T <C 1. The non-linear thermoelectric proper- 
ties of molecular junctions are therefore of great interest. 
Recent experiments [7|] probing the thermopower of thiol 
end-capped organic molecules showed non-linearities in 
the measured S already at AT w 0.1T. Earlier measure- 
ments of thermopower in metallic island single-electron 
transistors even displayed a change of the sign of the 
thermopower for very large AT jl5j ]. 

In the interesting regime of sharp electronic reso- 
nances, the electron tunnel coupling T is small and the 
main factor limiting the efficiency of molecular energy 
converters is expected to be the heat current from phonon 
exchange with rate 7 [10]. Nonetheless, to our knowl- 
edge, its effect has this far not been systematically in- 
vestigated. Only by making the tunnel coupling larger, 
T ^> 7, the phonon contribution to the heat current be- 
comes negligible. In this case, however, the efficiency 
becomes instead limited by the large electronic life-time 
broadening of the molecular resonances. The thermo- 
electric efficiency in this limit of coherent transport was 
studied very recently in the non-linear regime [16| us- 
ing both a many-body transport approach and a (non- 
interacting) approach based on Hiickel theory. Except 
for the latter work and a few others [TO, EJ , most theo- 
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retical studies of molecular thermoelectrics have focused 
on non-interacting models, using a Landauer type ap- 
proach. However, in the regime of weak tunnel coupling 
between molecule and electrodes, intra- molecular inter- 
actions typically constitute the largest energy scales of 
the problem. 

In this paper, we calculate the thermoelectric effi- 
ciency and converted electric power of a molecular device, 
including a single dominant molecular orbital, strong 
Coulomb interaction and coupling to a discrete vibra- 
tional mode, as well as coupling to lead phonons and 
lead electrons. Importantly, we include on equal foot- 
ing the phonon and electron contributions to the heat 
current, both of which contribute in establishing the sta- 
tionary occupation of the molecular vibrational mode. 
The Coulomb repulsion and electron- vibration coupling 
on the molecule are treated non-perturbatively in the 
limit of weak electron and phonon exchange in which 
thermoelectric efficiency is high. A central finding is that 
optimal thermoelectric operation typically is achieved in 
the non-linear, non- equilibrium regime. Here concepts 
of figure of merit and thermopower are no longer mean- 
ingful and the molecular occupancies, efficiency and out- 
put power must be explicitly calculated. The paper is 
organized as follows: Sect. [II] introduces the dissipative 
Anderson-Holstein model and the thermoelectric trans- 
port equations. In Sect, mil we present results for the 
efficiency and output power as function of the applied 
bias voltage and energy of the molecular orbital domi- 
nating transport. The heating of the molecule is ana- 
lyzed in Sect. HVl and the optimal choice of molecule and 
junction parameters is discussed in Sect. [Vj Section IVII 
summarizes and provides an outlook. 

Throughout the paper we set h = ku = e = 1, where 
h is Planck's constant, ks the Boltzmann constant and 
— e the electron charge. 



II. MODEL AND THERMOELECTRIC 
TRANSPORT THEORY 

Despite polarization and screening effects in molecular 
junctions [17), the electronic level-spacing in molecular 
devices is typically large compared to applied voltage- 
and temperature-bias. We therefore restrict our atten- 
tion to a single molecular orbital dominating transport. 
In fact, the thermoelectric properties have also been pre- 
dicted to be optimal in this case [l8|, [lj|. However, the 
quantized vibrational modes, which couple to the ch arg e 
localized on a molecular device, cannot be neglected |20l - 
[22) . The vibrations additionally couple to bulk phonon 
modes of the electrodes [23]. The goal of this paper 
is to clarify the importance of these excitations, char- 
acteristic of a molecular device, for the thermoelectric 
properties. We consider a thermoelectric junction as 
sketched in Fig. [TJ Its basic physics is captured by the 
following dissipative Anderson-Holstein model Hamilto- 



nian H tot = H + Hr + Ht , where 

H = en + Un^-ni + uutfb + uj\n (V + b) , (1) 
Hr = ^e kr c\ ar c k(Tr + ^2,u qr b\ r b qr , (2) 

kar qr 

H T = J2 trcl ar d a + h.c. 

kar 

+ Ecv^ + vM^ + fc)- (3) 

qr 

The molecular Hamiltonian, H, describes a spin- 
degenerate orbital level (operator d a for spin-projection 
a = ±|) with energy e and Coulomb repulsion U. The 
electron number n = n a , with n a = d\d a , is linearly 
coupled to the vibrational coordinate of the harmonic 
mode of frequency uo (operator b). The dimensionless 
electron-vibration coupling A is the shift of the vibra- 
tional potential as the molecule is charged, measured in 
units of the vibrational zero-point amplitude. The reser- 
voir Hamiltonian, Hr, describes the combined electron 
and phonon degrees of freedom in the two reservoirs, con- 
veniently referred to as the hot (r = h) and cold (r = c) 
electrodes. Non-interacting reservoir electrons with en- 
ergy e kr are created (annihilated) by c\ ar (c kar ); b\ r (b qr ) 
are the corresponding phonon operators for an electrode 
phonon mode with frequency uo qr . In each electrode elec- 
trons and phonons are assumed to be in equilibrium with 
temperatures Th = T + AT and T c = T respectively. 
The coupling between the reservoirs and the molecule is 
described by the Hamiltonian Ht, where the first term 
describes tunneling of electrons with amplitude t r and 
the second term couples the molecular and electrode vi- 
brational coordinates with (in general energy-dependent) 
amplitude C qr . In view of the thermoelectric efficiency 
we consider the case where both couplings are weak, i.e., 
we want a small tunnel broadening and a small heat 
current carried by the phonons. Therefore, Ht can be 
treated perturbatively below. The electron-electron in- 
teraction and electron- vibration coupling on the molecule 
are however allowed to take arbitrary values, which is a 
crucial aspect for addressing the important regime A ~ 1 
and U > T. 

We consider the thermoelectric junction in Fig. [TJa) 
operated as a heat to electric power converter. One elec- 
trode is heated (referred to as hot (h)) with the other 
electrode (referred to as cold (c)) kept at the ambient 
temperature. The hot electrode is grounded (chemical 
potential /i^ = measured relative to the electrode Fermi 
levels at zero bias) and a (negative) voltage —V is ap- 
plied to the cold electrode (/i c = V > 0). For simplicity 
we assume the capacitances associated with the tunnel 
junctions to both electrodes to be equal, resulting in a 
voltage-dependence of the molecular orbital, e oc V/2. 
We note that in an actual device which also makes use 
of the converted power, the voltage is not applied, but 
rather controlled by the temperature-bias and the resis- 
tance of the external circuit, see Fig. [U(b). 
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FIG. 1: (Color online), (a): Sketch of a thermoelectric junc- 
tion with a single molecule, drawn by way of illustration as 
Ceo- Operated as a thermal- to electric power converter, a 
temperature-bias is applied across the device. For an energy 
level above the electrochemical potentials, e > fi r , this can 
drive a net flow of electrons from the hot (h) to the cold (c) 
electrode by tunneling through the junctions. In addition, 
combined electron and phonon heat currents Qh and Q c are 
driven through the molecule, (b): The thermoelectric circuit 
is loaded by a resistor R, and as a result a voltage bias is ap- 
plied to the cold electrode, partially opposing the thermally 
induced electron flow. The voltage thus ranges from (corre- 
sponding to R = 0) to the value where 7 = (corresponding 
to R —> oo). In the rest of the paper we will however con- 
sider a test device as drawn in (a), where V, rather than R, is 
the free parameter. We investigate how to adjust V and the 
other parameters to obtain a maximal efficiency r/ and output 
power P. 

To formulate the transport equations, the linear cou- 
pling term in ffl is first eliminated by a standard trans- 
formation [24], H|| > which leads to a renormalization of the 
onsite and charging energies: e = e — A 2 cj, U = U — 2X 2 uo. 
After this transformation the eigenstates of H are easily 
found to be given by |a)|ra), where a = {0,t,|, 2} is the 
electronic state and m = {0, 1, 2, . . .} denotes the vibra- 
tional excitation number. The corresponding eigenener- 
gies are E am = E a + rnuo with Eq = 0, E a = e and 
E2 = 2e + U. Furthermore, the electron tunnel ampli- 



tude is renormalized to t r = t r exp[—X{b^ — b)], thereby in- 
corporating the Franck- Condon factors for electron tun- 
neling. The resulting transport characteristics under an 
applied voltage-bias have been analyzed in many works, 
see e.g., [25|-[27j. We note that in principle the additional 
coupling to reservoir phonon modes requires a more in- 
volved transformation, leading to more complicated ex- 
pressions for the renormalized parameters [23|. These 
corrections can be neglected in the regime considered 
here, where the coupling between the reservoir phonons 
and the molecular vibrational mode is weak. 

As mentioned above, the maximum efficiency of en- 
ergy conversion is expected in the limit of weak electron 
tunneling [lOj and weak coupling between molecular- and 
electrode vibrations: T ^> T r ,7 r . Here the rate for elec- 
tron tunneling involving electrode r is T r = 27rp r |? r | 2 , 
where p r is the density of states, which is assumed en- 
ergy independent (wide band limit). The relevant rate for 
phonon exchange with electrode r is 7 r = 27rd r (uj)\C qor | 2 , 
with qo defined by uo qor = uo. Here the phonon density 
of states, d r (E), and the coupling strength, C qri are in 
general energy-dependent. However, in the weak cou- 
pling limit only their value at uj qor = uo enters into the 
problem due to the selection rule \m — m'\ = 1 in low- 
est order perturbation theory in the coupling to elec- 
trode phonons (see the expression ([5]) for the rate ma- 
trix below). In the regime of non-linear temperature- 
and/or voltage-bias addressed in this paper, the molec- 
ular density matrix is not known a priori and needs to 
be calculated. For the weak coupling considered here 
this can be done using a standard master equation ap- 
proach. We can neglect contributions from non-diagonal 
elements of the density matrix, since the molecular states 
in our model are non-degenerate on the scale set by the 
rates (uo ^> T r ,7 r ) and spin-degeneracy does not lead to 
off-diagonal contributions. We note that this holds only 
in the weak coupling limit where the transport rates are 
evaluated to lowest order perturbation theory in T r and 
7 r [23. The transition rates for electron tunneling (W^) 
and phonon exchange (W^) can then be calculated from 
Fermi's golden rule [26l |28|: 
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^ r r |/ mm /| 2 {Sn a ,n a ,+l5M a -M^crfr(E arn — E a >m>) + 5 na }7la , _ \ $M' a —M a ,a [1 ~ fr(E a m — E a 'm')]} , (4) 
ar 

^7r<W {^m,m f +imb r (uo) + [&rM + l]}. (5) 

r 

I 



Here / mm > is a Franck-Condon factor [25], M a the spin- 
projection onto the z-axis of state a, which has electron 
number n a , and f r {E) = l/{exp[(E — fi r )/T r ] + 1} and 
b r (E) = l/{exp[E/T r ] — 1} are respectively the Fermi- 
and Bose distribution functions of lead r with electro- 
chemical potential fi r . The stationary state master equa- 



tion to be solved for the occupations P am then reads 

= E { W ala'm' + W^.J) P a >m>, (6) 

a'm' 

1 = ^P nm , (7) 



where Eq. ([7]) expresses probability normalization. At 
this point, we note that the interplay of charge and 
phonon tunneling is still non-trivial, as they do "inter- 
act" via the vibrational occupation number. A finite elec- 
tric current tends to highly excite the vibrational mode, 
leading to high effective molecular temperatures (see also 
Fig. |3(c)) and even clear deviations from an equilibrium 
(Boltzmann) shape of the distribution. This effect is 
particularly pronounced when the electron- vibration cou- 
pling is not too large, A < 1. The phonon current, on the 
other hand, tends to thermalize the vibration towards 
a temperature, that depends only on the temperatures 
of the hot (Th) and cold (T c ) lead and the relative size 
of the couplings 7^ and j c . However, through the exci- 
tations created by the electric current, the phonon cur- 
rent acquires an indirect dependence on both voltage and 
level position. Therefore accurate calculation of the non- 
equilibrium molecular state accounting for both electron 
and phonon effects is crucial. 

The electric current, 7 r , and heat current, Q ri going 
out of lead r, are given by 



Ed^kar 
dt 



= -EE(< e) ) , 

V / am, a'm' 
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dNko 
dt 



E< 



dM, 



J qr 



am a'm' 



EE «!) , , + « 

z — z — \ ^ / am, a'm' V 



where N^ar 



~ C kar C k< 



r and M qr = b\ r b qr . 



(8) 
(9) 

(10) 
P / / 

1 a'm' •> 

(ii) 

The electron 



dt 



am, am' 



current matrix Wj is similar to (J2J), but includes only 
processes involving reservoir r and a plus / minus sign 
for processes adding electrons to the molecules (first term 
m Q) / removing electrons from the molecules (second 
term in (J2])). Analogously, the heat current matrices, 
W^l and are similar to (|4]) and (j5]), respectively, 

but including only processes involving reservoir r and 
with the rate multiplied by the energy of the tunneling 
electron (measured relative to /i r ) or phonon. We note 
that there is in general some ambiguity associated with 
the definition of the heat current, see |29j, which however 
does not matter in the weak coupling limit discussed here. 

Beyond the linear regime, thermopower and figure of 
merit are no longer suitable quantities and we instead 
directly calculate the efficiency of the energy converter 
as follows. Driven by the thermal bias, electrons can 
gain potential energy by tunneling from the hot to the 
cold electrode via the molecule. The resulting electric 
output power is P = IV, where I = —Ih = I c . The 
input heat power, required to maintain the temperature 
bias, is equal to the heat current, Qh, flowing out of the 



hot electrode. The efficiency is thus given by 

V = IV/Qh. (12) 

Note that there is no conservation of the stationary heat 
current, as there is for the electric current, Ih + I c = 0. 
Instead the first law of thermodynamics guarantees that 

P = Qh + Qr. 



III. OPTIMAL BIAS VOLTAGE AND LEVEL 
POSITION 

We start by studying the efficiency and output power 
at fixed thermal bias, here chosen to be AT = T, as 
function of applied voltage-bias V and level position e. 
The efficiency of a single level quantum dot (spin-less 
electrons and no vibrational mode) was studied in [30| , 
where it was shown that the ideal Carnot efficiency is 
reached in the equilibrium limit of vanishing current, re- 
quiring the Fermi functions to be equal, /^(e) = / c (e), 
defining a line in the (V, e) plane: V = e(Th — T c )/T^. 
In Fig. [2] this equilibrium line corresponds to the bound- 
ary of the white areas. However, along this line also 
the output power vanishes (corresponding to reversible, 
infinitely slow operation without entropy loss). For van- 
ishing couplings to the phonon mode, 7 r — >• and A — >• 0, 
we recover this result in the non-interacting limit, U = 0, 
as well as for very strong interactions, U ^> T, AT. In the 
intermediate regime, the efficiency is slightly reduced. 

Switching on the electron- vibration coupling, but keep- 
ing 7 r = 0, the efficiency is decreased and never reaches 
the ideal value (ry c = 0.5 for AT = T), see Fig. C^a). 
In fact, 77 vanishes close to the zero electric current line 
(boundary of the white area), the reason being that, in 
contrast to the single- level discussed above, the heat cur- 
rent does not vanish completely when the charge current 
does. Inside the white area the current has been reversed 
by a too large voltage-bias and flows from high- to low- 
biased electrode and therefore does not accomplish any 
useful electric work (note that this regime can not be 
reached in the thermoelectric circuit of Fig. DJb)). The 
maximal efficiency is reached in the non-linear regime 
when the level is far above the Fermi edges of both leads. 
In this case electron transport involves very few ther- 
mally excited states in the heated electrode (tail of the 
Fermi function) and electron-induced vibrational excita- 
tions are exponentially suppressed, minimizing electronic 
heat loss. However, in this regime the current is highly 
suppressed, leading to a very small output power, see 
Fig. [2fd). Additionally, even a small coupling to the 
substrate phonons, 7 = T/10 in Fig. [2fb), drastically 
decreases the efficiency in this low-current regime, while 
having a much smaller effect in the regime where the 
current is larger (e is smaller). Thus, already a weak 
coupling to substrate phonon modes, 7 <C T, drastically 
changes the ideal operating conditions for maximum ef- 
ficiency by introducing a heat loss which depends only 
weakly on e and V (the dependence is indirect, through 
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FIG. 2: (Color online), (a)-(c): Efficiency 77 at thermal bias 
AT = T, as function of voltage-bias V and level position e for 
increasing coupling to substrate phonons, 7 = (a), 7 = T/10 
(b), 7 = T (c). In all plots A = 1, u = T, U = 10T and the 
couplings are symmetric, 7^ = 7 C = 7, I\ = T c = T. (d): 
Output power P as function of V and e for the parameters 
used in (b) (the power depends only weakly on 7) . 

the vibrational occupations). When the coupling to the 
substrate phonons becomes comparable to the tunnel 
coupling, 7 « T in Fig.(2fc), the efficiency is significantly 
decreased also in the high current regime. The output 
power, shown in Fig.[2fd) for the parameters used in (b), 
depends only weakly on 7 and is maximal for e « 4T and 
V « T. 



IV. TEMPERATURE DEPENDENCE AND 
MOLECULAR HEATING 

Next we fix the level position to a value with both 
large power and efficiency, e = 5T, and vary instead V 
and AT. The resulting efficiency and output power is 
shown in Fig. [3^a) and (b) respectively, for the same pa- 
rameters as in Fig. [2fb). As above, a too large voltage- 
bias compared to the temperature-bias reverses the cur- 
rent and no useful electric work is accomplished (white 
areas). The non-linear thermopower can be defined 
through V = -S(AT)AT at I = 0, i.e., V is the fi- 
nite voltage needed to compensate the temperature-bias 
and give zero electric current. Thus, —S(AT) is given by 
the slope of the line which passes through zero voltage 
at zero temperature-bias and hits the edge of the white 
areas at AT in Fig. 09(a) and (b). For large temperature- 
bias there are clear deviations from the linear response 
thermopower, S = S(AT — >• 0), given by the slope of the 
green dashed line in Fig. [3ja). 

As expected, the efficiency and (even more so) the 
power is increased by an increased temperature-bias. Fig- 
ure C9(a)-(b) shows the dependence on the temperature- 
bias over a wide range, all the way up to T^ = 3T. Such 
large temperature-bias could be obtained if the device 
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FIG. 3: (Color online). Efficiency 77 (a), output power P 
(b), and molecular temperature Tm (c), as function of V 
and AT, with e = 5T and other parameters as in Fig. [2Kb). 
Along the edge of the white areas in (a) and (b) I = 0, 
and this edge therefore defines the non-linear thermopower 
S(AT) = —V I AT at / = 0. The standard linear response 
thermopower, S = S(AT — > 0), is thus given by the slope 
close to zero, indicated by the green dashed line in (a). Note 
that Tm in (c) is plotted also in the regime where the current 
has been reversed, even though the device would normally 
not be operated under such conditions, (d): 77 (thick lines) 
and Tm (thin lines, normalized by 10T) as function of AT at 
V = 2.5T for 7^ = 7 C = T/10 (red solid lines), 107^ = 7 C = T 
(green dashed lines) , and 7^ = j c = T (blue dotted lines) . 



is operated at low temperatures. In applications, how- 
ever, T is most likely room temperature and junction 
stability limits operation to lower relative temperature- 
bias (e.g., AT = T/3 would mean T h w 400 K, which 
is a realistic value). However, it is actually possible 
to keep the (non-equilibrium) molecular temperature, 
Tm, much lower than the average electrode tempera- 
ture, Tm < T = T + AT/ 2, allowing operation at higher 
temperature-bias . 

To calculate Tm we use an idea suggested in [3l| and 
couple an additional phonon bath ("thermometer") very 
weakly to the molecule. The temperature of the ther- 
mometer bath is varied and Tm is defined as the bath 
temperature where the heat current between this bath 
and the molecule vanishes. Figure [3^c) shows Tm as 
function of V and AT, where it is seen that for small 
voltages Tm exceeds the average electrode temperature 
Tm > T. A larger voltage-bias, however, reduces the 
non-equilibrium electron current and Tm approaches T. 
Still, it is desirable to reduce Tm further, thereby al- 
lowing operation at higher AT without breaking the 
molecule. In designing molecular thermoelectric junc- 
tions it is therefore important to choose the electrode 
material and molecular anchoring groups such that the 
molecular vibration couples more strongly to the sub- 
strate phonons of the colder electrode (j c ^> 7^). This 
is shown in Fig. El(d), where j c = 7^ = T/10 (red 
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solid lines), j c = 10 = T (green dashed lines) and 
7c = lh — T (blue dotted lines). The asymmetric phonon 
coupling significantly reduces Tm (thin lines) by prevent- 
ing phonons from accumulating on the molecule: they 
enter slowly (7^) and exit quickly (7 C ). However, the 
asymmetry has a rather small effect on the efficiency 
(thick lines) since heat is still prevented from "leaking" 
through the molecule via the phonons as long as 7^ stays 
small compared to T. In contrast, with also 7^ large 
the efficiency goes down much more, and Tm goes up 
(the red solid and blue dotted lines for Tm still show 
a very small difference: the electron tunneling and sub- 
strate phonon couplings drive the system toward slightly 
different molecular temperatures, causing Tm to depend 
on the ratio 7/r). In realizing the mechanical coupling 
asymmetry j c ^7^, it is important keep the electronic 
tunnel coupling symmetric. Introducing an asymmetry 
in the electron tunnel couplings, Th ^ T c , while keep- 
ing + T c fixed, reduces the efficiency since the current 
level is set by the smallest coupling, while the phonon 
leakage current depends only weakly on and T c (only 
indirectly through the vibrational occupations). 



V. OPTIMAL THERMOELECTRIC 
JUNCTIONS 

Comparing Fig. [3ja) and (b), we reach an important 
result for optimizing molecular thermoelectric junctions, 
namely that for a given temperature-bias, maximum ef- 
ficiency and maximum output power is achieved at al- 
most the same voltage-bias. This is also seen in Fig. 2J 
where (a) shows the maximum efficiency, ?7max, and (b) 
the efficiency 77 at maximum power P = P max , both ob- 
tained by adjusting the voltage at given thermal bias AT. 
For non-zero coupling 7 to substrate phonons the max- 
imum efficiency and the efficiency at maximum output 
power are very close. The reason is seen from the rela- 
tion P = 77 (q£° +Q% ) ), where Q^ e) (Q^) is the electron 

(phonon) contribution to the heat current. Since 
only has a weak (indirect) dependence on the voltage- 
bias, 77 and P can be simultaneously maximized by ad- 
justing the bias when the phonon heat loss dominates 
(7 ^> T). As Fig. [3] and Fig. 2] shows, this holds approxi- 
mately also when 7 < T. 

When the electron- vibration coupling becomes strong, 
A > 1, the tunnel amplitudes involving the vibrational 
ground state become suppressed (Franck- Condon block- 
ade [lil). This reduces both the efficiency and output 
power since the current is decreased. Additionally, heat 
dissipation is increased since transport through excited 
vibrational states is favored, the typical energy trans- 
ferred to the vibrational mode by a tunneling electron be- 
ing given by the classical displacement energy, ojX 2 . This 
is shown in Fig. [5fa) and (b), where A = 2 (cf., Fig. [3ja) 
and (b)). Similarly, a smaller electron- vibration coupling 
enhances the efficiency, see Fig. Efc), where A = 0.5. 

In choosing uj = T 25 meV assuming room tern- 
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FIG. 4: (Color online), (a): Maximum efficiency, ?7max, as 
function of AT, i.e., the maximum obtained from vertical 
cuts in a plot such as Fig. [3]^ a). Parameters as in Fig. [3fa)— 
(c), but the strength of the coupling to substrate phonons is 
varied. Here 7 = 7^ = 7 C and T = Th = F c . (b): Same as 
(a), but instead showing efficiency at maximum output power 
i.e., 77 taken at the maximum P obtained from vertical 
cuts in a plot such as Fig. G^b). 




AT[T] 



FIG. 5: (Color online). Dependence on the molecular vibra- 
tion frequency uj and its coupling to the electron charge A, 
all other parameters fixed to those of Fig. EJa)-(c). (a)-(b): 
Efficiency 77 (a) and power P (b), as function of V and AT 
for A = 2 and uj — T. (c)-(d): Efficiency 77 as function of V 
and AT for A = 0.5 and uj — T (c) and A = 1 and uj = T/10 
(d). 



perature) we have investigated the influence of a rather 
high energy vibrational mode. Molecules often also have 
vibrational modes with much lower frequency, down to 
a few meV, especially when contacted to electrodes by 
linker wires [33[. However, as is shown in Fig.[5fd), where 
uj = T/10, a low-energy vibrational mode leads to a much 
smaller decrease of 77 compared to the ideal case of no vi- 
brational mode, as long as A is not too large. The reason 
is simply that a low frequency mode essentially can be 
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seen as a broadening of the electronic resonance of width 
~ cjA 2 , setting the scale for additional heat loss from 
electron tunneling compared to the case of no vibrational 
mode. Almost all decrease in efficiency in this case comes 
from the coupling to substrate phonon modes. In con- 
trast, a vibrational mode with a frequency much larger 
than the involved temperature- and voltage-bias, does 
not contribute at all to electron or heat transport (other 
than through the trivial shift of the electronic parameters 
e, U and t r through the electron- vibration coupling). 

Finally, we mention that in the simple model analyzed 
here, the strength of the Coulomb interaction does not 
play a crucial role. In the presence of a coupling to sub- 
strate phonons, reducing U leads to a somewhat larger 
efficiency (and output power) as the electric current is in- 
creased by the presence of another "transport channel". 

VI. CONCLUSIONS 

We have analyzed the efficiency and output power of 
a non-linear molecular thermoelectric device operated as 
a power converter. Accounting for the molecular vibra- 
tion and its coupling to substrate phonons turned out to 
be crucial in comparison with results for quantum dot 
models without these, as it qualitatively changes the op- 
erating conditions for optimal efficiency away from the 
equilibrium regime. By investigating a generic model sys- 
tem we can now identify some basic criteria for efficient 
molecular energy converters: 

(i) The coupling between substrate phonon modes and 
molecular vibrations should be asymmetric and mini- 
mal. Coupling more strongly to the colder lead re- 
duces the molecular temperature and allows operation at 
higher temperature-bias, improving efficiency and output 
power. 

(ii) The electron tunnel couplings should be symmet- 
ric, I\ & r c . Furthermore, they should be small as to 
minimize the life-time broadening, but still larger than 
the phonon coupling, T > 7. 

(Hi) The local electron- vibration coupling energy 
should be small compared to the zero-point energy of 
the vibrational mode (A < 1). This is most crucial for 
vibrational modes with frequencies around the operating 
temperatures and voltages. Modes with much higher fre- 
quencies do not contribute at all, and those with much 



lower frequencies only contribute to the heat loss through 
the coupling to substrate phonons. 

(iv) Ideal operating conditions for high efficiency and 
power is achieved when the conducting orbital energy 
is at a few /c#T from the Fermi edges of the electrodes 
(<~ 100 meV at room temperature). Control of the ther- 
mopower by adding electron donating or withdrawing 
groups to benzenedithiol molecules, thereby shifting the 
position of the HOMO and LUMO, was recently demon- 
strated @. The temperature-bias should be chosen as 
high as is allowed by molecular stability and the heat 
source. The ideal voltage-bias depends on the other pa- 
rameters, but is nearly the same when optimizing output 
power as when optimizing efficiency. Additionally, the ef- 
ficiency at maximum power is very close to the maximum 
efficiency. 

The general insights obtained in this exhaustive study 
of the most basic molecular thermoelectric model can 
serve as a guide for more complex molecular modeling, 
incorporating multiple vibrational modes, multiple elec- 
tronic states, breakdown of the Bqrn-Oppenheimer pic- 
ture (pseudo Jahn- Teller mixing (Hf), etc. In general, 
one expects deviations from a single-orbital model to give 
a less efficient energy-converter, as additional heat is lost 
by population of excited states. It might however be pos- 
sible to find special circumstances under which excited 
states can instead be desirable, e.g., by effectively cool- 
ing the vibrational mode. Atomistic studies of specific 
configurations of molecules, anchoring groups and elec- 
trodes may identify suitable systems which satisfy the 
above criteria and thereby further assist in advancing the 
chemical engineering of molecular thermoelectric junc- 
tions. For device applications, engineering of molecular 
monolayer devices, rather than ones based on a single 
molecule, presents a challenge to supramolecular chem- 
istry, nanodevice fabrication and surface science. 
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